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\0 ; Abstract 

O 

^^ ' We discuss the one-loop electroweak renormalization of the leptonic mixing matrix in 

the case of Majorana neutrinos, and establish its relationship with the renormalization 
group evolution of the dimension five operator responsible for the light Majorana neutrino 
Mh, masses. We compare our results in the effective theory with those in the full seesaw theory. 

^. ' 1 Introduction 

>■ 

p\ ' The complete renormalization of the Standard Model (SM) includes the renormalization of 

» . the mixing angles and phases of the CKM matrix. Given the experimental evidence for neu- 

trino oscillations, leptonic mixing occurs, and its pattern becomes part of the flavor puzzle. 
In the lepton sector, the renormalization group evolution of the flavor-mixing parameters 
may be quantitatively more significant than in the quark sector, as leptons mix with large in- 
tergenerational angles and the dominant gauge interactions are the electroweak interactions, 
which involve flavor mixing. In contrast, in the quark sector, the CKM matrix is close to 
diagonal and the flavor-blind strong interactions dominate. In this work, we discuss the rela- 
tionship between the one-loop electroweak renormalization of the leptonic mixing parameters 
and the renormalization group equations (RGE) for higher-dimension leptonic operators of 
the unbroken and spontaneously broken effective theories. 

Neutrino masses can be Dirac or Majorana. The Majorana case is natural in the technical 
sense of 't Hooft P, and is favored theoretically. Neutrinos are massless in the SM, since 
there are no right-handed neutrinos and the theory conserves {B — L), where B and L denote 
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baryon and lepton number, respectively. Non-null neutrino masses are a signal of physics 
beyond the SM, whose scale M is likely to be much higher than the electroweak scale. When 
M 3> Mz, the effects of new physics at the high energy scale can be parametrized in all 
generality by adding a tower of non-renormalizable effective operators to the SM Lagrangian 
which are SU{3) x SU{2) x C/(l) gauge invariant. The coefficients of the higher-dimension 
operators are suppressed by inverse powers of the large energy scale M. 

The lowest dimension non-renormalizable operator is the well-known dimension-five AL = 
—2 operator [5] responsible for neutrino masses, 



oif = Ul^0* uulp) , (1) 



where II^ ^-re the lepton weak doublets with fiavor denoted by the Greek index, the su- 
perscript c denotes charge conjugation^ and (p is related to the standard Higgs doublet by 
(f> = iT2<p* ■ The precise value of C5 is model-dependent. The leading candidate high-energy 
theory of neutrino mass is the seesaw model 0, the minimal extension of the SM which 
includes right-handed neutrino singlets Nr with heavy Majorana masses. The seesaw La- 
grangian is 

/:seesaw= = ill Ph + ie^PeR + iW^^Nr 



-h (t) Ye eR -iL^'Yv Nr - - Nr^ M Nr + h.c. , (2) 

where M denotes the Majorana mass matrix of the heavy neutrinos; and Ye and Y^, are the 
Yukawa coupling matrices of the lepton doublets to the right-handed charged leptons and 
neutrinos, respectively. In the seesaw model, the d = 5 coefficient at the seesaw scale is given 

by 

c5 = y;(M*)-iyJ. (3) 

Upon spontaneous breaking of the electroweak gauge symmetry (EWSB), this irrelevant d = 5 
operator yields a relevant d = 3 Majorana mass term for the left-handed weakly interacting 
neutrinos. The effective Majorana mass matrix of the light neutrinos is given by 

irnu)ap = -y (c5)„/3 (4) 

in the flavor basis. Diagonalization of the effective mass matrix results in flavor-changing 
couplings of W bosons to the leptonic charged currents. 

The RG evolution of the d = 5 operator coefficient C5 of the unbroken theory is given 
bym, 



■2 d 3 ■ 

167r — C5 = -- 

drnfj, 2 

+ 



C5{Y,yJ) + {Y,Y,Yc5 
3 g2 + 2 Ah + 2 Tr (3 Y^yJ + 3 Y^Y^^ + Y^yJ^'^ C5, (5) 



^The charge conjugate of the chiral fermion field is defined by ■tp'^ = 27072 V' • 



where Xh is defined as the couphng of the Higgs potential term —Xh {(!> (f) /4 in the SM 
Lagrangian. At order 1/M, the matrix C5 is the only source of lepton flavor mixing, and it can 
be diagonalized by a unitary matrix V, whose evolution has been studied ^. Since V is also 
the matrix which diagonalizes the light Majorana neutrino mass matrix in the spontaneously 
broken theory at order 1/M, Eq. Q), it can be identified with the lepton mixing matrix at 
this order. For instance, for the simple case of two light lepton generations in the absence 
of CP-violation, the lepton mixing matrix is parametrized by a single angle 6 whose RG 
evolution can be expressed in terms of the parameters of the broken theory. Eq. Q) yields 
the RGEs 
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where rrii and m2 denote the two eigenvalues of the light neutrino Majorana mass matrix 
Eq. (JH). 

The RGE of the lepton mixing matrix also has been obtained through the study of flavor- 
changing one-loop corrections to low-energy observables. Ref. [S] analyzes the case of the 
seesaw Lagrangian with both heavy and light Majorana neutrinos. In that case, the lep- 
tonic mixing matrix couples the charged leptons to both light and heavy neutrinos 7_. In 
consequence, the matrix which couples exclusively to the light neutrinos is not unitary. The 
one-loop counterterms for that matrix result in RGE with a rich structure which goes be- 
yond the approximations involved in obtaining Eq. @. For two light neutrino generations in 
the absence of CP violation, the RGE contains contributions not present in Eq. (jSJ above, 
proportional to 

2 2 
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as well as other terms present due to the non-unitarity of the mixing matrix. A term with 
mass dependence similar to Eq. ((7|) also appears in the RGE of the CKM quark mixing matrix 

m 

In this paper, we study the renormalization of the leptonic mixing matrix for the light 
Majorana neutrinos using the low-energy effective theory beyond the 0{\/M) approximation 
of Eq. ©. A term with the mass dependence shown in Eq. is quadratic in the light 
neutrino masses rrii and thus is an 0{1/M'^) effect from the point of view of the effective 
theory. A proper comparison with the full theory requires the derivation of the RGE in the 
effective theory to 0{1/M'^), which includes operators of dimension d = 6. The complete 
derivation of the one- loop RGE to 0{1/M'^) involves corrections from one-loop diagrams 
containing two insertions of the d = 5 operator or a single insertion of a d = 6 operator. In this 
work, we restrict ourselves to studying the dominant flavor-changing effects at 0{1/M'^). The 
Yukawa couplings Yi, are assumed to take natural values of 0(1) throughout. All corrections 
at order 0{1/M) are included, but at order 0{1/A'P) corrections suppressed by Y^ and 
gauge couplings are neglected. We also determine the relationship of these results to the 
counterterms obtained in the full theory in which both light and heavy neutrino species are 
present. 



One could ask why a similar analysis is not necessary for the neutrino oscillation proba- 
bility, and more concretely for its Am^ dependence, which is an 0{1/M'^) effect. The answer 
is that the dominant one-loop corrections to the neutrino oscillation probability are those 
affecting the neutrino propagator, and while terms quadratic in C5 give contributions to it 
in the seesaw model, terms linear in the d = 6 operator coefficients do not. In contrast, 
the definition at one loop of other low-energy observables, such as the W decay width into 
leptons, does require inclusion of all d = 6 operators, and their one-loop mixing with terms 
quadratic in the d = 5 operator or linear in the d = 6 operators. In this sense, such low- 
energy observables depend on additional aspects of the new physics at high energy, which are 
encoded in the d = 6 operator coefficients of the low-energy effective theory. 

In section 2, we obtain the RGEs for the effective theory before electroweak symmetry 
breaking (EWSB): the evolution of the d = 5 operator is derived in the R^ gauge to 0{1/M), 
and the RG evolution of the dominant flavor-changing contributions to 0{1/M'^) is deter- 
mined. Section 3 performs the same analysis after EWSB. Section 4 compares the results 
for the effective theory with the counterterms obtained in the full theory in which heavy and 
light neutrino fields are present. Feynman rules and details of the computations are included 
in the Appendices. 

2 RGE before EWSB at 0{1/M^) 

The low-energy effective Lagrangian of the minimal seesaw model to 0(1/M^) is of the form 

Ces = CsM + S£''=^ + 6£''=^ + ■■■, (8) 

where £sM is the SM Lagrangian and the higher-dimensional operators give the low-energy 
physics effects of the heavy Majorana neutrinos. 

At the seesaw scale, the d = 5 term of the effective theory is given by Eqs. (^ and (jH)). 
The d = 6 term is [H] 



^aP 



Oif ^ [iLa<t>) ip [j^^^Lp) , (9) 

where 

c6 = y.(|M|2)-iyJ. (10) 

When the Higgs doublet acquires a vacuum expectation value, this irrelevant d = 6 operator 
leads to a relevant d = 4 flavor-nondiagonal kinetic energy term for the left-handed Majorana 
neutrinos. After normalization of the kinetic energy term and diagonalization of the kinetic 
energy and mass terms, the d = 6 operator ultimately results in a non-unitary correction to 
the leptonic charged currents and light Majorana neutrino neutral current 0. 

At the seesaw scale, the coefficients of the operators O and O are given by Eqs. © 
and ()1U|). Additional d = 6 operators will be present in the low-energy Lagrangian, however, 
since they are generated by radiative mixing in the renormalization group running of the 
d = 6 effective Lagrangian from the high-energy scale M to the low-energy scale //: 

5£'^=6(/x) = 6C^=\fi)^ = ceifi) 0'^=6 + J2 [c^(/^) Oiifi) + Q(/i)t 0,(/i)t] , (11) 



where i runs over all independent d = 6 operators. 

For our purposes, we only need to consider those operators which are not flavor blind 
and contribute to the lepton masses/kinetic energies and to the leptonic mixing matrix upon 
EWSB, that is, to those operators which ultimately modify either the charged and neutral 
currents or the charged lepton Yukawa coupling matrix Y^, up to 0{1/M'^). 

First, consider the one-loop corrections to the d = 5 operator coefficient, the 0{1/M) 
term. The complete one-loop renormalization of C5 has been studied in the literature in the 
t'Hooft-Feynman gauge 1^. We have re-derived the results in the covariant Rg gauge. The 
bare and renormalized O5 coefficients are related by 
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where 2e = 4 — d, and Z^ and Z^^ are the one-loop wave-function renormalization constants 
of the scalar and lepton doublets in the MS minimal subtraction scheme: 
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From Eqs. H12 |) - (|15|) . one obtains the RGE of the d = 5 operator given in Eq. Q. 

Now, consider 0{1/M'^) corrections. The d = 6 operators do not correct O 
0{1/M'^), as there are no SU{3) x SU{2) x U{1) gauge invariant and (B — L)-non invari- 
ant d = Q operators [H]. The coefficients of d = 6 operators are corrected at 0{1/M'^) by 
effects quadratic in the d = 5 operator and by effects linear in the d = 6 operators. We 
will restrict our analysis to 0{v'^ /M"^) and OiY^ v^ /M'^) corrections, disregarding in this 
work 0{g'^ v^ /JVP), 0{g''^ u^/M^) and 0{Y^ v^ /M"^) contributions. There are no corrections 
proportional to the Higgs self-coupling \h to 0{v'^ /M'^). 

We classify the possible d = 6 operators by using the operator basis defined in Ref. pTTj . 
with the following exception: the weak singlet and triplet operators 



are replaced by the orthogonal combinations 
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d=6 
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(18) 



This replacement is advantageous because 0\^ is equal to the d = 6 operator O'^ ^, Eq. ©. 

The one-loop diagrams resulting in relevant 0{1/M'^) corrections from two insertions of 
the O operator are given in Fig. 1. These diagrams result in contributions to several 
operators defined in the basis of Ref. ^3 • Details of the computation of the diagrams in the 
Feynman gauge are given in Appendix ^ The contributions relevant to the flavor-mixing 



renormalization program modify the coefficients c 



(-) 



ce and c\« corresponding to the 



operators in Eq. ([T8|l. and the coefficient Ce</, corresponding to the operator 



O. 
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defined in Ref. PTO]. As will be discussed in the next section, O^,^ is important for our analysis 
since it results in fiavor non-diagonal corrections to the charged lepton masses, and thus to 
flavor mixing in charged currents when the Higgs acquires a vacuum expectation value. OL ' 
also is kept as it contributes to flavor-changing charged and neutral currents. 

The one- loop diagrams (see Appendix^ result in the counterterms: 
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where the counterterms are deflned to remove the divergent 1/e poles of the one-loop graphs 
in the MS scheme. 
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Figure 1: One- loop diagrams with two insertions of the d = 5 operator. The extra arrows next 
to lepton indices indicate fermion flow, as deflned in Ref. 1111 



The bare and renormalized coefficients of the operators OL = O'^ ^ , OL and Oec/) are 
related by 

4?'^^° = (4)"'' «)"'' Vt^ + ^4t^] /^" (^^.)-^/^ (^.)"^^ ' (23) 
c^r = {^y'" ^z^y'" {^l)"^ [^^^ + ^^-t>\ /^'^ (^<^)"'^' (^e«)-'/' , (24) 

Eqs. ^-^ lead finally to the RGEs: 
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These are the results before EWSB. It is necessary to consider the spontaneously broken 
theory in order to compare with low-energy experiments. 



3 RGE after EWSB at 0{l/M'^) 

There is an important difference between the effective mixing matrix defined at tree-level and 
at one- loop, when EWSB is taken into account. At tree- level to 0{1/M'^), all charged and 
neutral-current mixing effects stem from flavor non-diagonal terms in the neutrino masses and 
neutrino kinetic energies, which arise from the operator coefficients C5 and cg of the effective 
theory. At one-loop to ©(l/M^), however, additional effects arise from flavor non-diagonal 
contributions to the charged lepton masses. 

3.1 Tree-level effective Lagrangian 

After spontaneous electroweak symmetry breaking, the relevant part of the d < 6 effective 
Lagrangian is given by 

>Ceff = i eia P GLa + i eRa P GRa - {cZa {jne)al3 GRp + h.C.) 

+ iVLaP{5al3 + Kp) l^LfS " -^ {j^a {■my)al3 l^Lf3 + h.C.) (28) 



+ -^W.Z + -^His 



;s ) 



where 



Cw,z = -^ {J-''''W+^ + J+^'^'W-^) + -^J^^Z^ (29) 

Y 2 COS c/p^ 

contains the W and Z couplings to the weak charged and neutral currents, and >CHiggs contains 
the lepton couplings to the Higgs boson. 



In the flavor basis, the charged and neutral currents are given by 

J'^^^ = eLo'lfi'^La, (30) 

The charged lepton mass matrix, the light Major ana neutrino mass matrix and the neu- 
trino kinetic energy flavor-changing effective coupling at tree-level are given by 

rriv = — ^C5, (31) 

A = yc^. =yC6. 

As shown in Ref . ^9^ , the effect of the d = 6 operator coefficient A is transferred to the 
weak currents by redefining the neutrino field to the canonical basis after EWSB (or by using 
the equations of motion of II before EWSB). The net effect is a rescaling of the neutrino field 
i^La ^ (I — 2)013^^113 in the effective Lagrangian Eq. (|28|) . After the rescaling, the charged 
and neutral currents are given by 

cc _-_ .. ( X _ {^Cc)cx(3 



Jf!^ = l^T^Lal^ {Saf3 - {Kic)ap) ^Lf3 + ( "o ~ '^^^'^ ) ^^"^A^^ia + {-^'^^w) eRal^idRa , 

where 

^cc = AJVc = ^ (33) 

at tree-level. We will see in the following subsection that Xcc and AJ^^^ have different one-loop 
corrections. 

It is customary to define the tree-level matrices in the flavor basis where nie is diagonal, 
but rrii, and A are not. The mass basis is obtained from the flavor basis by making a unitary 
transformation V which diagonalizes the light Majorana neutrino mass matrix, 

V*m,V^ = {m,)^,,^^^^,,. (34) 

The light neutrino Majorana mass eigenstates Vi = vf, i = 1, • • • , n, are defined by 

ly^ = V,aULa + V*,UL''^, (35) 

SO 

l^La = Pl (v^) .Vi = PL VL n . (36) 



ai 



The weak currents in the mass eigenstate basis 

J-^"""" = eZ^l^.Kf'^i, (37) 



are defined in terms of the non-unitary effective leptonic mixing matrices 

ub'h = f I - ^) y^ , (38) 



<c - (l - ^) Ft . (39) 

At tree-level, Uf^ = ^NC 

Recall that rriu does not get a contribution from A up to order 0{1/M^). (The rephasing 
of the charged lepton fields in Ref. |2j is ignored here.). 

3.2 One-loop effective Lagrangian 

The one-loop effective Lagrangian is given by 

£^r = ^eff + SCs (40) 

where C^s is the renormalized Lagrangian and SC^s is the counterterm Lagrangian which 
cancels the ultraviolet (UV) divergences of the one-loop diagrams. In the canonical flavor 
basis, the counterterms relevant for one-loop leptonic mixing are 

SCeS = i {5Ze^)^p eLapeLjs + i {5Zej^)^p eRapejifs - (ela (5me)Q,/3 e/j/3 + h.c.) 

+ ^ i^^i^L )a/3 T^La P T^Lf3 " ^ (^a {^^u)al3 ^Lp + h.C.) (41) 

+ -j= [cLa 7a. (5Acc)a/3'^L/3 W^^" + h.C.) 



COS Ow 



-I^Lalu i^^Nc)aP ^Lp " ^^^aT/i {^^Nc)aP ^Lp 



Coupling constant counterterms bg^ 56 w are not relevant for the mixing analysis and have 
been omitted. 



In the MS scheme, the bare and renormalized quantities are related by 

^Bare ^ (^ J-1/2 [^^ + 5^^] ^2e ^-1/2 ^ (43) 

m^-^ = (ZjJ-'/' [m, + 6m,] fi'^ ^-^2 ^ (43) 

>^cF = {Ziy'^' [Xcc + 6Xcc] Z-ll-" , (44) 

{y^c?""'^ = [ziy" [A5Vc + <JA5^c] zzl'^ , (45) 

(A^^)B- = (4^) "'^' \\%c + ^^%c\ Z-l'^ ■ (46) 

In general, the one-loop wavefunction renormalization constants Zg^ and Z^j^ are modified 
from Zi^ in Eq. H13|l by the diagrams in Fig. 2: 

Ze, = Z,^-r^\M'^C%\ (47) 

Z., = Z,,-r^\{Ml + Ml)c^;,\ (48) 
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Figure 2: One-loop diagrams with one insertion of the d = 6 operator. 

In our analysis, we neglect these additional contributions since they are 0{g'^v'^ /M'^). For 
completeness, we also give here the wavefunction renormalization constant of the right-handed 
charged lepton in the R^ gauge: 

The vertex counterternis are given by 



5mu = -— 5c5 — C5 , (51) 



(52) 
(53) 
(54) 





2 " 2 " ' 


5^cc = 


v^ ^i'i^ - ^c^^i^ 


2 2 


^K^c = 


2 0« ' 


^K^c = 


2 '^'^'/'^ • 



Notice that flavor-changing corrections to the left-handed electron neutral current are 
generated at one-loop. 

3.3 RGE of U^% 

The operator coefficient Jcg,/, induces flavor-nondiagonal one-loop corrections to the charged 
lepton mass matrix, so it is necessary to diagonalize the charged lepton mass matrix at one- 
loop. The unitary matrix E which diagonalizes the charged lepton mass matrix is defined 
by 

E (meml) E'f = (m^)diag,real • (55) 

The relation of E to its bare matrix is 

^Bare = f = E + 6E . (56) 
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The presence of a non-trivial matrix E changes the definition of the charged-current lepton 
mixing matrix, 

J;:'''' = eIalAUccUi^i, (57) 

which is defined at one-loop by the renormalized non-unitary leptonic mixing matrix 

Ace' 



Ufc^ElI 



yt 



(58) 



Note that the product EV' is the leptonic analogue of the CKM matrix. 
The RGE of U^ is then given by 



""^-nU^l-U^m-lE^V^ 



a in ^ 2 a m ^ 

where Tg and Ti, are anti-Hermitian auxiliary matrices introduced to maintain the unitarity 
of E and V under renormalization group running, 

dE 



dln/i 

dV 

dlnfi 



TeE, 



T^V. 



Te can be derived from Eq. (|55|) . 

V^^e/diag,real ^ -'e i^^e Jdiag,real \^^e) 
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Te+E 
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dln^ 
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with the RGE for the charged lepton mass matrix given by 
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(60) 

E\ 
(61) 

■ (62) 



Imposing the constraint that the off-diagonal matrix elements on the RHS of Eq. (J6T 
vanish to order 0(1/M^), it follows that 



{Te)a(i 



, a = /3 



(63) 



where the non-diagonal terms stem from the RGE of the operator coefficient Ce<^. 
A similar expression can be found for T^, |l2j : 



Re [{T,)ij] 
Im [{T^)ij] 





1 3 g^ m^+nij „ 
16^ 2 2M^ nii-mj ^^ 



1 3 g^ m.-nij J 
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I = J 
i = j 



(64) 
(65) 
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Finally, the last term in Eq. (|59|) can be rewritten as 
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using Eqs. ^, ^^ and (fSU)) . 
Thus, Eq. ((^ becomes 



dln/x 
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For illustrative purposes, let us write explicitly the result for two generations of light 
fermions in the absence of CP violation. The effective mixing matrix at tree level is then 
parametrized by 



u: 



eflf 
CC 



1 _ Ace -^EM 

^2 2 



2 



'MM 

2 



cos a sm f 
— sin 9 cos I 



(68) 



The corresponding RGE of the matrix elements {U'^q) . are 
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3.4 RGE of Ufc 



From the definition of Uf(j in Eq. (jSHJ, we obtain its RGE at ©(l/M^) in terms of the 
running of the unitary matrix V in Eq. H60() and the running of \^^(j'- 
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where Tj^ is given by Eq. (|E5|) and 
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Thus, Eq. (|71|1 becomes 
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m^. + m^ ) Vjp nip V*p + 2mjmi V*^ rrip Vip 



(73) 



For the example of two light lepton generations considered previously, the RGE are 



dlnfi 



9 



32tt'^M^ 2 






(m2 



2 2 ^""2 



mo — m^ 



ml) 



sm f cost 



H cos t' TTi]^ > , 



d Wc) 



e2 



(iln/i 



327r2M2, 2 



?TT,^ + TTT-g 



mg 



mf,-m^ 



(mi + 7712) 



2 v^a ~ "T-i) "I ^2 _ ^2 



TWo — TTlt 



(m^ - m^) 



smf cos 



sm t/ rrin 



(74) 



4 RGE in the High-Energy Seesaw Theory 



Consider now the full seesaw Lagrangian in Eq. Q. The results in Ref. obtained with n' 
right-handed neutrino generations and n lepton generations can be compared with the RGE 
obtained in the low-energy effective theory at 0{1/M'^). We begin by reviewing the high 
energy theory. More details can be found in Appendix C. 

After spontaneous breaking of the electroweak gauge symmetry, the portion of the seesaw 
Lagrangian (jJJ which breaks the chiral symmetries of the lepton kinetic energy terms is given 
by 

1 



rxSB 

■^seesaw 



-^ulY^Nr + — ei Ye BR + h.c. 



{Nr'^ MNr + NrM*Nr') .(75) 



The left-handed neutrino fields can be arranged in a (n + n') column vector 



riL 



Nr' 



Eq. (|75|) can be written in terms of ul as 
rxSB ^ 



1 V 

- ul" MnL + —T^eiYe cr 



+ h.c. , 



(76) 



(77) 



where the [n + n') x [n + n') neutrino mass matrix M is defined in terms of the nxn' Yukawa 
coupling matrix Yi, and the n' x n' Major ana mass matrix M of the heavy neutrinos. 



M 



JLV* 

% Y^ M* 



(78) 



V2 
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The neutrino mass matrix M is diagonalized by a {n + n') x [n + n') unitary transformation 
V, 

V*A^Vt = (A^),i,g^_i, (79) 

where the diagonal mass matrix is described by its {n+n') mass eigenvalues. Light eigenstates 
Ui get mass eigenvalues nii, i = l,---,n and heavy eigenstates Nj get masses M/, / = 
n + 1,- ■ ■ ,n + n' . The n + n' Majorana mass eigenstate neutrinos n = n^ are defined by 

^^[n) =VnL + V*<. (80) 

It is convenient to define submatrices of A^diag, real and the diagonalizing matrix V |13j : 

v.(--). ,S, 

The neutrino weak eigenstates are 

i^La = Pl ( [k^) ^^ n + (X^) ^^ iV/) , (83) 

so that the weak currents are given by 

^;'''' = ^7M(i^i^^)a^^*+^7M(^^^U^^/, (84) 

J^'~^ = ( -- - S^Qw ] eLaliJ.eLa + {-S^Ow) eRalfieRa (85) 

+ \ (17, 7^ (i^ K^)^^ V, + Nj 7^ (X X^)jj Nj + 17, 7^ (i^ Xt)^^ iV^ + Nj 7^ (X i^t)^^. ^^.) , 

where E' is the unitary matrix which diagonalizes the charged lepton mass matrix mg. 

As our interest is in the renormalization of the weak currents which couple to the light 
neutrino eigenstates, we define the non-unitary matrices 

Ucc = EK^ (86) 

Unc = K^ (87) 

whose RGE will be compared below with the corresponding matrices obtained earlier in the 
effective theory at 0(1/M^). To this level of approximation, one can solve for the submatrices 
of V and A4diag,reai from the constraints that V be unitary and that it diagonalizes the mass 
matrix M. The complete analysis is given in AppendixO The solution to 0{1/M'^) for the 
mixing matrices is 

K = v(l- ^xA = ^ ^ - ^ 
X = -V*hx\ 
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where the matrix x [H] is defined by 



X 



V2 



Y^I 



-1 



(90) 



The matrix V in Eq. (|88|) is the same n x n unitary matrix which diagonahzes the hght 
Major ana neutrino mass matrix at 0{1/M'^), previously defined in Eq. (|34() . 



V* (m^) . 



(91) 



'^diag, real 

The matrix Vh is the n' x n' unitary matrix which diagonahzes the complex, symmetric, 
heavy Major ana neutrino mass matrix at 0{1/M'^), 



Mm 



diaK,real = Vh [ M* + ]^(f X* M* + M*xh 



V, 



H 



(92) 



Thus, to ©(l/M^), 

Ucc = EK^ = U§j , (93) 

Unc = K^ = Ufc , (94) 

with U^(j and Uf^(j being the mixing matrices of the effective theory, Eqs. (|38|) and (|39)) . 

The RGE of Ucc and Ujyc have been computed in Ref. .0! in the on-shell (OS) renor- 
malization analysis of the CKM matrix [S] . In order to compare the results in the full seesaw 
theory "W with the results in the effective theory, we expand the former to 0(1/M^) using 
Eqs. H88|) and (|89() . In the basis where me and M are real and diagonal, the RGE of Ucc is 
given by 

,2 



d{Ucc)c 
dln/Lt 



dlnn 

+ e{M-fi) 



1 1 g' 



V* m'^ 
167r2 2 2M^ '" ' 



e{pi-M) 



167r2 



m^ + "^a 



2 ■^-^ rrin 



mS 



(95) 



{YuY})^p iUcc)pi 



3 v:^ 



{Ucc)aj 



2 ^ m^ 
and the RGE of t/Arc is 



2 



(m^ + mj) Vjp (y.yJ)/3^ 1/4 + 2m,-m, V*. (r.^J);^ F^ 



7/37 ''n 



d([/. 



NCJai 



dliifj. 



dlnfi 



e{fi - M) 



1 7 5" , .* 2 

— V ■ Tfi- 

167r2 4 2M,2 *^*° '"^ 



'VK 



+ e{M-fi) 



1 



3 -^ 

+ 9E 



167r2 

{UNc)aj 



-Kp{YuYl)p^,V, 



J7 



(96) 



jV« 



TTT-^ 



2 

raj 



im] + m?) V,-^ (i;yJ)/37 1/4 + 2m,-mi K*. [Y,Y%^ V 



7/37 "n 



The 6 functions in the equations have been introduced hy hand to single out those terms which 
are different in both theories. The differences between the two theories stem from diagrams 
with heavy neutrino eigenstates running in the loop. Depending on the effects below the 
scale M at which the heavy neutrino eigenstates decouple, we can divide these diagrams into 
two classes: 
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Diagrams which are UV divergent in the full theory and zero in the unbroken effective 
theory. For instance, the term 



3 " m|+m| 
present in Eq. (|95|) . stems from the diagram in the full theory in Fig. 3. 



(97) 



Full Theory 


Effective Theory 


/ \ 

/ \ 

► 1 1 - 


/' *'" '\ 

\ / 

\ / 

\ / 

* ^-O-^^ — * 

eL/3 CLoL 


eL/3(p) Ni 


'eLa(p) 


IGtt^ e (^^^^ )"/5 2 


ig"^ 1 X „ *^ 



Figure 3: Matching self-energy diagrams in the full and effective theories 

The corresponding diagram in the effective theory, where the heavy neutrino propa- 
gator is replaced by a non-local interaction given by the O"^^^ operator, vanishes in 
dimensional regularization in the unbroken theory. 

In a mass-dependent regularization scheme, the diagram is not zero; instead it is propor- 
tional to the square of the cutoff which cancels the 1/M^ factor of the d = 6 operator. 
This apparent contradiction between mass-dependent and mass-independent treatments 
is solved when all orders of perturbation theory are resummed in the mass-dependent 
scheme, recovering the null result obtained in dimensional regularization. After EWSB, 
this diagram does contribute positively to the RGE evolution, however, as the figure 
indicates. 

All contributions arising from diagrams of this type are proportional to Y^Yj . Another 
example is provided by the term 



1 



-A„^(y,yJ)^^F4 



(98) 



in Eqs. ()95() and H96|). which also vanishes in the unbroken effective theory, although it 
contributes after EWSB at 0{l/M^), see Fig. 4. 



16 



Full Theory 


Effective Theory 






Viiv) Nj Vk ; Ni 

1 
% 


1 


i^i(p)' " z^fc /"\ 'vi{p) 
X X 


ji^\{VY,Yiv%i 


0(M-3) 



Figure 4: Matching self-energy diagrams in the full and effective theories 



Diagrams which are UV finite in the complete theory and diverge in the effective theory. 
In addition to the terms in Y,^Yj explained above, the terms 



167r2 2M^ 



Vni mf 



(99) 



have a different numerical weight in Eqs. (|95jl and (|96jl . The origin of this difference 
is the same as for the case above: the replacement of the heavy neutrino propagators 
by the non-local interactions below the decoupling scale. Here, the heavy neutrino 
propagator is inside a non-UV divergent loop and its removal increases the degree of 
divergence, yielding UV divergent terms which contribute to the RGE. An example is 
given by the diagram in Fig. 5 showing a contribution which is UV safe in the full 
theory, but which is divergent in the effective theory. The same occurs in Fig. 6 where 
only one of the four diagrams which contributes to the neutrino self-energy diverges 
in the full theory: in the effective theory a factor 4 appears as the four diagrams are 
all UV divergent after decoupling. By studying the possible combinations of external 
legs, one realizes that the first three diagrams contribute to the effective operator OL , 
which also contains couplings V ecj)'^ (fP* ., see Fig. 6. The last diagram in Fig. 7, which is 
UV-safe in the full theory, gives rise to a divergent coefficient for the effective operator 



O 



(1) 



Full Theory 


Effective Theory 


/ \ 


/ \ 

ep{p) 1 \ eiaip) 

— ^ ^- 

1 1 
1 1 
k k 


Nj \ Vk \Nj 

1 1 
1 1 


UVsafe 


it^ 1 i^c^u i 



Figure 5: Matching self-energy diasrams in the full and effective theories. 



Full Theory 



Effective Theory 



f 



^ -N 



Vj (P) 



Nj Vk Nj 



i^iip) 



i 1 



T6^ e 



y{y^i,yJ)(y:i^y^)y' 



i 

2 



U,{p) ,Nj , Nj Ui{p) 



k 



k 



UVsafe 



Uj (p) Nj I Z^fc Nj ^ Vi (p) 



I 



UVsafe 



Uj (p) Nj ^ Vk ^ Nj Vi (p) 



X X 



J^j [p) I 

I 
I 

k 



Vk 



I Vi (p) 



UVsafe 



4T^7(c5C5)a/3 



Figure 6: Matching self-energy diagrams in the full and effective theories 
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Full Theory 



Effective Theory 



T^j (P) 



eaip) 



Nj 



^k 



Nj 



I 



-i 1 
T6^ e 



{y^^y.l {y^ik^y^) y\ 



aj_ 



■C'uv 



l^j (p) I 

* O- 

I 



\ 



T^k 



^3 (P) 



eiaip) 



\ eLa{p) 

-^ 

I 

I 



Nj Uk Nj 



I 



I 



UVsafe 



-^i2(c*5C5Wf 



Figure 7: Matching self-energy diagrams in the full and effective theories 



5 Conclusions 

We have studied the renormalization of the leptonic mixing matrix in the case of light Ma- 
jorana neutrinos. The analysis has been performed in the low-energy effective theory to 
0{1/M'^), with M being the generic high-energy scale at which lepton number is broken. 
The leading flavor-changing RGE have been obtained including one-loop corrections which 
involve two insertions of the d = 5 neutrino mass operator as well as tree-level contributions 
of d = 6 effective operators. 

In the existing derivations of the RGE at 0{1/M), the matrix which diagonalizes the 
d = 5 operator coefficient in the unbroken theory also diagonalizes the neutrino mass matrix, 
and thus is the physical mixing matrix (when working in the basis in which the charged 
lepton Yukawa couplings are diagonal, as is customary). In contrast, one-loop effects at 
0{1/M'^) corresponding to two insertions of the d = 5 operator coefficient result in flavor- 
changing corrections to the charged-lepton mass matrix. Its diagonalization also affects the 
evolution of the physical leptonic mixing matrix, just as in the quark case. This realization 
is crucial to the comparison of the effective theory results with the renormalization program 
of high-energy theories in which both heavy and light neutrino fields are active. 

Furthermore, the tree-level and one-loop insertions of the d = 6 operator characteristic 
of the seesaw theory yields non-unitary contributions to the mixing matrices of the effective 
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theory, for both charged leptonic and neutrino neutral currents. 

We have finally compared our results in the effective theory with the RGE of the mixing 
matrices of the full seesaw theory developed to 0{1/M'^). The terms which differ in both 
expansions have been understood consistently in terms of the different degree of divergence 
of the diagrams in both theories. 
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A Feynman rules for d = 5 and d = 6 operators 



The d = 5 effective Lagrangian for the seesaw model is given by 



= ^{C5)ap {hl4>*) {4>^hp) + 1^{CIU (^J) [^^(^l'p) 

where a,b,m,n are SU{2)l gauge indices, and e = iT2 is the totally antisymmetric tensor in 
two dimensions. The Feynman rule for the d = 5 vertex is: 




The extra arrows next to lepton lines indicate fermion flow, as defined in Ref. jll| 



e{k) ei-k) ^ 

The factor of 1/2 in the Lagrangian term is cancelled by combinatoric factor of 2 for 
a <-> /3 of the two lepton lines. The two factors of ee gives the combinatoric factor of 2 for 
two (j) lines. The coefficient C5 satisfies c^ = C5. The vertex is AL = —2. 

The Hermitian conjugate d = 5 vertex is: 
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i {cDad (e""" e^ + e'^" e^) (^) 



'^Laa 



The vertex is AL = 2. Notice that chirality projection operator is now Pr. 
The d = 6 effective Lagrangian for the seesaw model is given by 



5C' 



d=6 






The Feynman rule of the d = 6 vertex is: 



(100) 



4>n{p-q') \ / (t>mip-q) 



^Lp,{q' 




i{c,)^pe---e^^{i>) (i^) 



^Laaiq) 



The heavy right-handed neutrino which was integrated out had momentum p''. The cg 
coefficient satisfies cl = cq. The vertex is AL = 0. 



B |c5p Diagrams 

The BW d = 6 operators induced by one-loop diagrams with two insertions of the d 
operator are: 



O 



^^} = iihal^hp) {(t>^D''<p) 



o 



(3) 



i{iLa7^.r£Lp) (ck^D^^rcl) 



O 



(1) 



o 



o. 



(3) 
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OjJ^ = ^ {iLal^.iL^) (iL^l'^iLs) 



and the Hermitian conjugate operators. 
We will use the operator set: 



of = (0t(D^<^))((D/^0t)^ 
Oe4. = {iL<peR) (^<PU) , 



so that Oq is an operator in the basis. 

The effective Lagrangian is defined to be 



C,s = Y,CiO,, (101) 



which implies that the operator Oi produces vertex with Feynman rule iciS, where S is the 
symmetry factor of the operator. (Equivalently, the operator Oi/ S, defined with its proper 
symmetry factor, produces vertex with Feynman rule icj.) 

Coefficients cL = cg and c,^ are related to cU and cL by 

'-(fit — ^(j)E ^(j)E ' 
J+) _ ^(1) 1^(3) 

and 

'^(t>e ~ 2 \ vi ^ <t>^ 
^(3) - \(S+)_S-) 

The UV parts of the one-loop diagrams with two insertions of the d = 5 operator are 
calculated below. 
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B.l 



</>« {p - q') 



hp,{q') 




4>m (p - q) 



La 



a{q) 



"Te;;^! 2 (c5*C5)a/3(i') y 2 
X ( A f<^'"^ fbn _|_ !:abfimn\ 



The first term in parenthesis yields a contribution to the coefficient of Og 

1 1 



Acs 



2 (C5*C5) . 



(102) 



167r2 e 
The second term in parenthesis gives rise to the following combination of d = 6 operators: 



i {ha 7m hp) ( 0^ D^'ct^ - {D''<P^) <P) + i (ha P hp - ha P h^ ( ^V) 



.(103) 



The first term is the operator OJ with a coefficient given by 



(1) 



Ac 



(1) 



and it gives contributions to 



Acit) = AciV 



16^ 7 4 ("^ ''^ 



1 1 1 



Ac6 = Ac' 



(1) 



167r2 e 4 
1 1 1 



- T IC5 C5j , 



167r2 e 4 
Using the equations of motion of the lepton doublet 



(C5*C5) . 



i p lL-Ye(t>eR + ■■■ = Q, 

-ill p -e^<A^yJ + --- = o, 



(104) 

(105) 
(106) 



(107) 
(108) 



where the ellipsis stands for terms suppressed in 1/M, the second term generates Oecf, with a 
coefficient 



Ac, 



it-W^'"^''''^^ ■ 



(109) 
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B.2 



'"Laa 



4>n{q) 




/ <t>k{q') 



^ I xnl Tmk I x'K'k ^iml^ 



16 



^ \ Tr[c5C5*] {p + qY {d^'^S"'^ + 5'^^8' 



i>m{p) / 
/ 



"-LPh 



(t^iiv' =v + q- q') 



There is an explicit symmetry factor S = ^ associated to this diagram. The diagram 
contributes to the following combination of operators 

h\A /"^t n2^^ _L (a,] 2)2^^ ^^]^ +2i (t)^ D^'A { (j)^ Df'ct) 1 



D' 



+ 



2 (^ (^t J:,^^A f (DM0t) (p\+2((l)^A( (L»^(/.t) D^'A 



+ 



D' 



b^D^ 



The first term on the right hand side of this equation is 20\ . The second one is 20\ . The 
third and fourth terms on the right hand side can be rewritten using the equation of motion 
for (p, 

D^^ = d^YjqL + qZ^yuUR + e]^Y^UL (HO) 

and shown to cancel. 

Operators O, and O, result in the Feynman rule 2 (^,5"'^™^ + ^nk^mi'-j ^ggause of a 
symmetry factor S = 4. Considering all symmetry factors, the following coefficients are 
induced by the above diagram 



Ac' 



(1) 



1 1 



Ac' 



(3) 



167r2 e 
1 1 



167r2 



Tr[c5C5*] 
Tr[c5C5*] 



(111) 
(112) 



B.3 




'"Laa 



LPb 



levr 



^ 7 (c5)a/3 (£5)7^ [^Lpl^tLs] [^Lalf^ih^) 



Again, there is a symmetry factor S = i^ and a symmetry factor S" = 4 in the Feynman rule 
for the operator [iipl'^^Ld) {^LaliM^L'y) ■ Thus, the above diagram generates the coefficient 
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since 



is defined with 1/2. 



Ac 



O 



a/B'jS 



_ 1 

yjS'yS 2 



Y^72^^^^"^ ^^^^^^ 



{harhp) (iL^l^ils) 



(113) 
(114) 



B.4 




i 1 



ml rafc i rmfc ral \ 
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^i (c5*C5ye)a;3('^™'<5'^'+'^'"'^'5' 



^Laa 



eR/3 



The above diagram contains two identical external (p hnes, so there is a symmetry factor 
5 = 2 in deducing its coefficient from its Feynman rule. Thus, the above diagram induces 
the operator Oes with a coefficient 



ACe0 = -T-^ - {c5*C5Ye) 



(115) 



Notice the difference with the case of diagram Bl. There, we obtained Oe^ from applying 
the equations of motion to another operator which does not contain identical particles. Thus, 
we didn't apply a symmetry factor correction in Eq. (7). 

B.5 

In summary, the one-loop diagrams with two insertions of C5 give a divergent contribution to 
a number of the operator coefficients. Specifically, 



AC6 



^^t' = -l^K' + i'^^'^^' 



Ac 



Acp 



(+) 



1 1 1 



167r2 e 4 
1 1 3 



C5 C5. 



C5*C5 Ye] 



Ac 



(1) 



167r2 e 4 

1^ 7 ^'^''*''^ ' 



(116) 

(117) 
(118) 
(119) 
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^^6^ = TT^-Tr[c5*C5], (120) 



167r2 



e 



,(,)\ 1 1 1 



or 

^'''^^ " 167r2 e "^^ ""' 

C High Energy Seesaw Theory 

Considering Eq. I\8'2\\ . unitarity of V implies 

KK^ + WW^ = I, (122) 

XX^ + ZZ^ = I, (123) 

XK'^ + ZW'' = 0, (124) 

K^K + X^X = I, (125) 

W^W + Z^Z = I, (126) 

W^K + Z^X = 0, (127) 

so that the submatrices K, W, X and Z are not unitary. 
Eq. ((7^ imphes that 

mdiag.real = -^ [k*Y:W^ + W*Y^ K^) + W* M*W\ (128) 

= ^ (k*Y*Z'< + W*Y^X'<] + H^*^^*^^ (129) 



Mdiag.real = ^ ( X^Y,* Z^ + Z*YJ X^) + Z* M* Z\ (130) 



V2 

One can solve for the submatrices of V and A^diag.real from the constraints that V be 
unitary and that it diagonalizes the mass matrix Ai. The solution to 0{1/M'^) is 

K = v(l-^xx^^=v(l-^y (131) 

W = Vx, (132) 

X = -V*hx\ (133) 

Z = V*h{i-\x^^. (134) 

as weh as Eqs. (J91l92p . 

The matrix V to 0{1/M'^) is given by 

y=(V{l-\xX^) Vx \ , . 
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